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Abstract 



We consider a homogeneous gas of spin-S" fermionic atoms, as might occur 
near the center of an optical trap. In the case where all scattering lengths 

> 

C ' are negative and of the same magnitude we demonstrate the instability of 



the Fermi sea to the condensation of bound 'baryonic' composites containing 



O 

Q\ ' 25 + 1 atoms. The gap in the excitation spectrum is calculated. 

£ ■ I. INTRODUCTION 

C 

Sympathetic cooling of trapped fermionic gases M is likely to lead to the creation of 

O ■ 

degenerate Fermi gases in the near future. In magnetic traps, the fermion's spin is locked to 
>< : the field direction. However the advent §g of optical traps fields the promise of degenerate 

Fermi gases with unquenched spins and hence a variety of nonideal Fermi gases. Already 
there has been some discussion of the possibility of pairing in the presence of attractive 
interactions, which may be of an exotic type |7|||. In part the exotic possibilities occur 
because alkali fermions exist which have larger total spins than spin 1/2: 22 Na, 86 Rb, 132 Cs, 
134 Cs, 136 Cs have S = 5/2, 5/2, 3/2, 7/2, 9/2 respectively §. In this paper we will point 
out that in addition to pairing there are other possible ground states, in particular when 
the atoms experience spin-independent attractive interactions. 

One-dimensional models of interacting fermions with spins greater than one half have 
been considered for a substantial time, both for repulsive [0] and attractive [pd|-[T3|] inter- 



1 



actions. In the attractive case, the ground state has been found to contain bound states of 
fermions with up to 2S+1 constituents (with equality in the absence of an applied 'magnetic' 
field which would distinguish the different spin states). One may rationalise these results by 
noting that the Pauli principle does not militate against binding extra fermions until a spin 
state must be doubly occupied, and hence extra nodes in the spatial wave function occur. In 
the light of these results we will examine whether the ground state of the t/iree-dimensional 
homogeneous S > 1/2 weakly interacting attractive Fermi gas is characterised by a conden- 
sation of composites of 2S + 1 fermions. The possible condensation of alpha particles as 
four-particle composites in nuclear matter and at the surface of nuclei has a long history 
whose relation to the current work we will discuss at the end of this paper. 



II. BARYONIC GROUND STATE ENERGY 

In three dimensions we use a variational approach similar to that used in the original 
BCS paper |TJj]. We will assume that the composites ('baryons') have zero centre of mass 
momentum (as in BCS) and are also total spin singlets, as in the one-dimensional results. 
(The latter can be understood physically as providing the lowest kinetic energy associated 
with relative motion in the bound state.) The condensed state, with all baryons having 
centre of mass momentum zero, is of the form (here n = 2S + 1): 

N 



IB) 



^2 <f(h, ■ ■ ■ , k n )cl ltai ■ ■ ■ cl njan 



where all the spin states in the baryon, <7j, (with i = 1, • • ■ , n) are different. ip(ki, • • ■ , k n ) 
is the Fourier transform of the completely symmetric relative wavefunction (p(ri, • " ,r n ). 
Unlike the BCS case, there are several grand canonical states, |BCS), corresponding to |B). 
We will construct the simplest one. 

This family of variational states allows us to use 1/n as a small parameter. (Indeed if 
the coupling between the atoms were strong, there would be similarities with work |T5] on 
the structure of baryons in QCD using a 1/N expansion.) 



Consider the following (non-normalized) ground state: 



\ip h ) = exp 



y{k u ...,k n )c\(ki)...c ] n {k n ) 



10} 



Below we show that the corresponding normalized state has the form 



|V>b) = e 



(n-2)7V b /2 TT 



Y\u n ~ l (k) - exp 



(fa(h, ...,&„) Cj(/Ci) 



10} 



>-2)jV b /2 



n 

fc! + ... + fc n =0 



+ . . . , fe n )ci(A;i) . . . 4(fc r , 



1=1 



10} 



(1) 



where c^(fc), c a (/c), a = 1, . . . , n = 2S+1 the creation and annihilation operators of a fermion 
with a-th projection of the spin and momentum k, N\> = N/n is the number of baryons 
(with iV being the number of atoms) and 



u(A;i, . . . ,k n )u(k a ) 

(p a {K\, . . . , K n J = 7~» \ n x ffci+...+fe„,0 



(2) 



w(fei) • . . . • u(fc n ) 

The Kronecker symbol on the r.h.s. of @) means that we consider baryonic states with total 
momentum equal to zero. The expression k\ + . . . + k n — as a subscript to the product 
in (|l]) implies that the product is taken over all sets {k\, . . . ,k n } with the total momentum 
equal to zero. We will apply a 'normalization' condition similar to that used in BCS states: 



u 



(£;)| 2 +K£;)| 2 = 1 , \v(k)\ 2 = £ \v(k, k 2 , . . . , K 

k2,—,k n 



(3) 



This, unlike the BCS one, does not imply that the total state is normalised. The vector 
conjugate to is 



(^ b | = e (»-%/2J| u n-l W . (Q| exp 



<fa(ki,...,k n )c n (k n )...ci(ki 



(4) 



To calculate the normalization (V'bl'^b) we make use of the following identity 



exp 



Y Paikl, ■ ■ ■ , knjclikt) . . . 4(/c n ) 



fci .....kn 



D£D£ exp 



i=l k 



(5) 



where are Grassmann variables and the measure of integration is denoted by 

n-l 



8=1 fc 



As enters linearly in the exponentials in (|5|) then the operator averaging in (ipb\ipb) can 
be easily fulfilled and we find 



<Vb|Vb> = C/ DCD^ exp 



n-l 



i=l A; 



2 <fa(h, • • • , k n )ip a (p 1} . . ■ ,p n )Cl(Pl)6(^l) • • • Cn-lCPn-OCn-lC^n-l)^^ 



where C denotes the factor 



For the second term in the exponential we make use of the following identity: 



exp 



- 5Z <Pn(fa, ■ ■ ■ ) ^nVn(Pl, • • • ,Pn)fl(Pl)fl(fcl) • • • £n-l(Pn-l)£n-l(&n-l)4„,P« 
Pl,--,Pn 

DaDcr exp 



n— 1 ra— 1 

53 5Z^( fc >p) ff i( fc 'P) - 53Z^^(^'^)^(^)^(^) 

i=l fc,p i=l k,p 

n-l 

+ 5Z <Pn{fa, ■ ■ ■ , k n )<Pn(Pl, ■ ■ ■ , Pn) Yl a i(Pi' k i)$k n ,p 



1=1 



(6) 



Pl, -;Pn 



where 



n-l 



DcrDa = JJ JJdcr i (A;,p)dcr i (A;,p) 

i=i fc,p 

The meaning of a-fields is not clear at this point, but we will discuss this presently. 
Integrating over the Grassmann fields we obtain 

n— 1 n— 1 

53 53^ (k, p)ai {k, p) + 53 trlogfl + Gi] + K 

i=l k,p 



i=l 



(V'blV'b) = C / DaDcr exp 
where .fT denotes the following expression 

K = 5Z ^n(fcl,---,^n)^n(Pl,---,Pn) II a i{Pii ki)5k n ,p, 



n-l 



8=1 
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We calculate the integral over a, a using the saddle point method, where n is the large 
parameter. We assume that the symmetry of the exponential over <7j-fields is not broken 
and put (jj = a, Oi = a, so the exponent becomes 



(n - l)aa + (n - l)tr log[l + a] + K 



The saddle point equations are 



a(k,p) = (l + a)-\k,p) 



n— 1 



a(k,p) = ( Pn{k,k 2 ,...,k n )(p n (p,p 2 ,...,p n )Y[cr(pi,k i )6, 



(7) 
(8) 



fc 2l ...,fc n 

P2>--->i>n 



i=2 



In the limit n — > oo we expect a mean field approximation to be valid for the description of 
the composites: 



v(h, ...,k n ) = B- 5 kl+ ... +knfi J[ v(ki 



(9) 



i=l 



where B is chosen to satisfy the normalization condition (|3|). With this assumed form for v, 
it can be shown that the following anzatz can serve as a quite general solution of the saddle 
point equations: 



a(k,p) = S^pS^k) + c 2 s 2 (k)s 2 (p) , cr(k,p) = 5 k , p s 3 (k) + c i s i (k)s i (p) . 



(10) 



with functions Si,i = 1,2,3,4 to be determined. If we further restrict ourselves by consid- 
ering weak coupling case we obtain that the solution is (see Appendix for details) 

v(k) 2 



a(k,p) = 5 k)P \u(k)\ , a(k } p) = 5 k , p 



u(k) 



(11) 



Hence, 



in 



-l)J2Hk)\ 2 -(n-l)^og\u(k)\ 2 + Y: E \v(k,k 2 ,...,k ri 

fc A> /c ^2 5***1 ^ ?i 



(il>h\ipb) = Cexp 

Because of the normalization condition @ and expression for C we finally obtain 



(V>b|Vb) = l + 0(l/n) 



(12) 



Let us consider the following average: 

(A\c a (k)ci(k')\tfj h ) 

Repeating all the steps we have made before we can obtain the following result: 
(A\c*(k)cl(k')\^ h ) = C J DaDa (1 + a a )~ x {k, k') 



x exp 



n— 1 n— 1 

i=l fc,p i=l 



Calculating again the integral over a-fields using the saddle point method and noting that 
the contribution of the first multiplier in the integrand to the saddle point equations can be 
neglected we obtain from (0): 

(^\c a (k)ci(k')\^) = a(k,k') 

So the cr-field may be interpreted as momentum distribution of the holes. 
In an analogous way we can obtain the following averages: 

(i; h \ci(k)c a (k)\i; h ) = \v(k)\ 2 + 0(l/n) (13) 

E(^b|4(A;)4(-A; + q)c p {-k' + q)c a (k')\^ h ) 
i 

= <WK&)| 2 E N<?)| 2 + v(k)u(k)u(k')v(k')f(k, k') + 0(l/n) (14) 
g 

where k = -Mj and the function f has the form 

\k\ ■> 

/(fc,fc') = /(x) = l-y + y , X = Sm 6 - (15) 

with 9 being the angle between vectors k, k' . 

We are now in a position to consider a general Hamiltonian of the form: 

H = J2< k )cUk)c a (k) + ±- J2 V{k,ky a {k)c\{-k + q )cp{-k' + q)c a {k') (16) 



2fi m*„ 



where 



<*)-£-/. (17) 

and is the chemical potential. Then, using the results above for expectation values, we 
find the following expression for the ground state energy 

k "' " k,k' 

(18) 

with the condition 

(^ h \N\tfj h ) = n^2Hk)\ 2 = N (19) 
k 

Note that if we consider an interaction potential depending only on the modulus of mo- 
mentum: V(k,k') = V(\k\,\k'\), and look for a solution for u, v also depending only on 
the modulus of momentum (that is in absence of a spontaneous breaking of the rotational 
symmetry) then the function / in (|I8]) can be replaced by its average value: 

i 3 
f(x) = - 

and we have 

E h = (Vfc|#|Vfc> = nJ2e(k)\v(k)\ 2 + + ^ £ V(k, k')v{k)u{k)u{k')v{k>) (20) 

where V = |V. 

The normalization condition, (|^), allows the introduction of the following parameteriza- 
tion for u, v, by analogy with the usual procedure for the BCS case: 

v(k) = cos9(k) , u(k) = sm6(k) (21) 

and we obtain for the ground state energy 

V(0)N 2 n 2 

E h = nY, e(k) cos 2 9{k) + y ' + — £ y (fc, k ') cos d ( k ) sin #(^) cos 9(k') sin 6{k') (22) 

, AH ZS I ■ , / 

We now minimize with respect to to find 



tan[20(fc)] 



1 n 



2e{k) (2tt) 3 
Introducing the following notation 



dk' V(k,k') sin[20(fc') 



(23) 



n 



2{2nf 



dk' V(k, k') sin[20(fc / )] , E(k) = e 2 (k) + A 2 (k) 



1/2 



(24) 



such that 



tan[20(Jfe)] = , sm[29(k)\ = ^® , cos[20(fc)] 



e(k) ' E(fc) 
we obtain the following equation for the gap A(k): 



E(k) 



(25) 



A(k) 



n 



2(2tt) 3 

and the chemical potential can be defined from 



dk' V(k, k 



E(k') 



(26) 



7? 



(2tt): 



/" dfc' cos 2 6{k) = p 



(27) 



Following Anderson and Morel | lql , Eq.(p4|) can be rewritten in terms of an effective 
potential, U^(k,k'): 



A(k) 



1 



/ dk'U £ (k,k 

J\e(k')\<£ 



(2vr) 3 J\e{k')\<Z 



n A(A/) 

2£(A; , ) 



(28) 



where £ is some cut-off such that £ <C /i = ^ and t/g(/c, fc ) satisfies the equation 



£^(fc, fc') = nV r (ife,A; / ) - 



1 



dq nV{k,q)U c {q, k') 



1 



(2vr) 3 J Hq) \>t " y 5V " '2|e(g)| 
We take K(/c, fe') to have the following separable (energy dependent) form: 

Ana 



(29) 



V(k, k') = VQ(n - \e(k)\)e(n - |e(fc')|) , V 



rn 



(30) 



which, as we shall see, is consistent with low-energy approximation to the T-matrix. Here 
a is the scattering length (we put H = 1), 0(x) the step function. We assume that the 
interaction is weak (which is a good approximation experimentally) which means 



k-pa -C 1 
8 



(31) 



To check that the form of interaction (p0|) is a consistent low-energy approximation consider 
the equation for the T-matrix T(k, k', z) of zero energy, z = : 

i /* 777 

T(k, k', 0) = V(k, k') - -—3 / dq V(k, q)T(q, k', 0)- (32) 
Making the anzatz 

r(fc, fc', 0) = T 6(/i - |e(A;)|)©0 - |e(fc')|) 
we obtain the following relation for T , V: 

1 1 m£; F 

p T p V ' P ° 2tt 2 

As 

1 7T 

> 1 



2akp 
then 

T « V (33) 
Let us define U^(k, k'). Assuming the form 

U^k,k')=U^(p-\e(k)\)e(p-\e(k')\) 



we get the following equation for U^. 



nVUf f q 



2 



4tt 2 J * |e(g)| 

$<|6( 3 )|< M 

After some algebra we arrive at the relation 

-L = — ^ + log(^2 + 1) - log i (34) 
PoUt; PonV 4/i 

If we assume again that we can drop the second term in the right hand side (because the 
coupling is weak) and noting that V, < we finally obtain 

' 1 ^ + log 7- (35) 



Po\ u t\ Pon\V\ 4/i 
Note that the energy can be expressed in terms of 



E h = — — - / dA; 

2(2tt) 



e(k)- 



e 2 (k) A 2 (k) 
E(k) ~ 2E(k) 



(36) 



From (28) then we obtain 

A = 2£exp 

The energy of the normal state per unit volume can be written as 



1 


= 8/i exp 


1 


Po 






p n\V\_ 



(37) 



E n 



71 



dk e(k) 



(38) 



Hence, from (^f 
energies is: 

AE h 



(27r) 3 J\k\<k F 

]) we find that the difference of the baryonic and the normal ground state 



Eb — E n 



-48p/i exp 



P 



X 



Po 



mki 



(39) 



p n\V\_ 

with iV the number of atoms in the system. 

We conclude this section with a note about weakening the assumption of equality of 
scattering lengths in all channels. If the scattering lengths are different then the interaction 
will take the form | jT7|j : 



[/] 

V = 2 V n (S 1 ■ S 2 f 

n=0 

where [/] means the integer part of /. In the case where the scattering lengths in all channels 
are equal, a F = a, WF, we have V = —, V n = 0,n ^ 0. If the scattering lengths are slightly 
different then the terms in the interaction with V n ,n ^ are small and can be neglected 
but Vo is equal to some weighted average over the scattering lengths and the treatment in 
this section will still be approximately valid. 



III. BCS GROUND STATE ENERGY 

We will now show that the conventional BCS ground state has an energy that is higher 
than the baryonic one, at least if we restrict ourselves to s-wave pairing. The BCS ground 



state has the form [18 



10 



I^bcs) = II u n (k)exp 

k>0 



10} 



(40) 



and the conjugate state is 



(v>bcs| = n^ n ( fc )(°i ex p 



fe>0 



- v{k)c a {k)P a pc p {-k) 



k>0 

a, f3 



(41) 



where 



v(k) 
u(k) 



\v{k)\ 2 + \u{k)f 



(42) 



k > is an arbitrary ordering on momentum space which divides it into two halves (it can be 
defined, for example, as follows: k > if k z > or k z — 0, k y > or k z — 0, k y — 0, k x > 0). 
The matrix P has the form 



-P 



(43) 



T means matrix transposition. 

It is well known that [I5| any antisymmetric matrix A : A T = —A can be orthogonally 
transformed to an antisymmetric matrix B in the canonical form: 



A = OBO 1 



o T = o- 1 



matrix B has the following quasi-diagonal (block-diagonal) form 

Ai OA 



B = diag ^ 

-Ai — A n / 2 



n/2 



This transformation is equivalent to 



c(Jfe) -> Oc(k) , c\k) -> c^ife)© 1 



(44) 



As the Hamiltonian in (^) is invariant under the transformation ( f4"4| ) and because of the 
assumption of the spin symmetry of the ground state, there exists some orthogonal matrix 



O such that the ground state (j40| , |4TD , after transformation (|44]) , takes the following form: 
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I^bcs) = II u n (k)exp 

k>0 



E <p(k)4.(k)Paf*i(-k) 



10} 



(45) 



where 



P = dia£ 



e 1Xl 



_ e ixi o ' -e ix "/ 2 



(46) 



The form of the matrix P implies that the spin states can be enumerated after the trans- 
formation in such a way that there are n/2 "positively" directed spins and each of them 
has an oppositely directed partner. So the spin indices can be thought as taking values 
±1, . . . , in/2. The phases Xa represent possible phase differences between Cooper pairs in 
different spin states (we consider only rotationally invariant states, so the modulus of all 
elements in the matrix P are equal). 

To calculate the norm and the energy of the ground state we consider the following 
generating functional: 



G(rj, -q) = (?/>Bcs|exp 
Making use of the following identity 



k.a 



cxp 



k,a 



I^BCs) 



(47) 



cxp 



k>0 

a,0 



J D£D£ exp 



E L(k)Uk) 



k>0,a 



+ E ci(k)Uk) + J2Uk)P a M-k) 

k>0,a 



k>0 
a. 



we get the following result for the generating functional 



G(n,n) = exp 



n/2 



E E {m<k)jX"f) a (k)fj„ a (-k) - u(k)v(k)e- iXa V«(k)V- a (-k)} 



k a=l 



-EEI#)l 2 ^(%( fc ) 

k a 



(48) 



Putting n,n = in (|48| ) we obtain (^bcsI^bcs) = 1- The energy of the ground state is 



^BCS = (V'BCSI^I^BCS) 



n E e(k)\v(k)\ 2 + + ^E V(k, k')u{k)v{k)u{k')v{k') 

k ' " k,k' 



(49) 
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Introducing as usual the trigonometric parameterization 

u(k) = sin9(k) , v(k) = cos 6{k) 

and applying the variational principle we find the following equation 

1 1 



tan[20(fc)] 



2e(fc) (2tt) 3 



dk' V(k,k') sin[20(fc')] 



(50) 



Defining the gap as 



A(k) = -- 



dJfe' V(k, k') sin[20(A/) 



(51) 



2(2vr) 3 

and repeating calculations of the preceding section we arrive at the following expression for 
the gap 

1 



A = 8/zexp 



Po\V\ 



(52) 



and the difference between BCS and the normal ground state energies is 



AE) 



BCS 



E 



BCS 



n 

From eqs.( 



n 



—A8pfi exp 



1 



Po\V\ 



P 



X 

n 



Pa 



mki 
~2~^ 



(53) 



5^) we conclude that the energy of the baryonic ground state is lower than 
for BCS state (for n > 4 or, equivalently, S > 3/2). Indeed, one usually assumes that a 
coupling constant is proportional to 1/n in the framework of large n expansion, so we can 
put V = Vq/u. Then the ratio of the baryonic and BCS binding energies can be expressed 
as follows: 

AEbcs 



AE h 



exp 



1 1 

TFH + ^7 

Pol V^j Po^l^l 





3n/8 - 1" 


= exp 


Pol Vo| 





< 1 for n > 4 



The last relation demonstrates the statement. 



IV. CONCLUSION 



In this paper we have shown that, assuming that all scattering lengths are approximately 
equal, the true ground state of a dilute gas of fermions with a high hyperfine spin is in fact 

13 



of baryonic nature. The energy of the state and the gap in the spectrum (interpreting A by 
analogy with the BCS case) have been calculated. For comparison we calculated the energy 
of BCS s-wave ground states and showed that it is higher than the energy of the baryonic 
ground state. 

The possibility of forming a baryon-like bound state has been discussed in paper |EU 



The authors investigated the formation of a four-fermion (a-particle) condensate. While the 
results of |20| are very interesting they cannot be directly related to ours. Firstly, in the case 
of a-particles n = 4 and 1/n = 1/4 can be hardly regarded as a good expansion parameter 
to apply our results. Secondly, n = 4 is the dimension of the joint spin-isospin space and 
the Hamiltonian is not invariant under rotations in that space (exemplified by the absence 
of a di- neutron bound state, as against the existence of the deuteron). 

We note that if the interaction were strong then the quasi-baryons considered in the 
paper would become 'real-space' composite bosons where the ground state might be more 
appropriately described as being Bose-condensed. So it would interesting to study the 
evolution from weak to strong coupling in the manner that Nozieres et. al. derived for an 
attractive fermion gas |2T| . 



The method developed in this paper may be of interest in other fields such as the low- 



energy behavior of QCD |22| , |23| , the quark-gluon plasma |24J] and neutron stars [|25H27 
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APPENDIX: 



In this appendix we sketch the main steps of calculation of the norm of the baryonic 
ground state vector and averages ( |T3|JT^ ). 

We consider the following generating functional 
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G{r],rf) = (ip h \ exp 



£E4(*H(*) 

.i=i k 



exp 



i=l k 



l^b) 



(Al) 



from which we can obtain the norm and the required averages in an obvious way. Using the 
identity (13) we can transform (|A1|) to the form 



G(f}, V ) = C f D£D£ exp £ + Vi(k%(k) + ^k)ru(k) - f}i(k) Vi (k)} 

J U=i k 

+ X {&n(h, ■ ■ ■ ,k n )r] n {k n )^i{k n -i) . . -Ci(^i) 
+ip n (k 1: . . . , fc„)fi(fci) • • •^n-l(^„-l)^„(A; rl )} 



(A2) 



- X! ¥>»(fcl, • • • , U^n(pi, • • • ,Pn)fl(pi)fl(fcl) • • -Cn-lCPn-OCn-lC^n-l)^,^ 

Here we put a = n while it can be any number from 1 to n. 

Consider first the norm of the ground state. Then after some transformations (see the 
main text) we obtain saddle point equations (|7|||). If we accept anzatz (0) then from (||) 
we have 



s 3 (k) = K 



v(k) 



u(k) 



s 4 (/c) 



v(k) 
u(k) 



n—l 



K = K Q (ki) = \B\ 2 £ 5 kl+ ... +knfi H 

k2,-.,k„ 



i=2 



vjki^s^ki 
\u(ki)\* 



\v{k n )f 



- (I \ I D|2 r x TT V ( k i) V (Pi 

C 4 = C 4 ^l,PlJ — \&\ 2^ (5 fe 1 +...+fe n ,Ofpi+-+P»-l+fer l ,0 11 x / 

P2>--->Pn-l 

rn-i 



X 



n-1 



IH^P^^) + c 2S 2 (Pi)s2(^)} - II $k l)Pl Sl{ki 



(A3) 



.i=2 i=2 

K (k),c^(k,p) are weakly dependent on their arguments and can be approximated by con- 
stants. On other hand, from (0,0) we obtain 



si(k) 



l + s 3 (k) 



c 2 = -c 4 



s 2 (k) 



s 4 (fc)P 



Choosing I? such that K = 1 we obtain the following solution 

\v(k)\ 2 



(A4) 



Sl (k) = \u(k)\\ s 3 (k) 
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\u(k)\' 



S2(k) = v(k)u(k) , s^k) 



v{k) 



It can be easily seen then that the condition Kq = 1 is equivalent to the normalization 
condition From (|A3| , [A4|) a relation for determining of the constant C4 can be obtained. 
In the weak coupling approximation it takes the form: 



const 
n 3 N h 



k 



;i + c 4 iv b )™- 2 

An obvious solution of this equation is C4 = and we obtain (IT)- 

Consider now average (p!3|). Using the symmetry over spin we obtain from (S3) after 
integrating over the Grassmann fields: 



(il) h \c\(k)ci(k)\i) h ) = {i) h \c ] n (k)c n (k)\tp h ) = - 



d 







dr) n (k) drj n (k) 



G{r),rj) 



n-l 



ri=-q=0 



C J Do-Da J2 <Pn(h,---,k n -i,k)tp n {pi,...,p n -i,k) Y[[l + CTi] l (k i} p. 

fci.— .*n_i 7=1 



x exp 



PI.— 

n— 1 ri—l 

^2^2a i {k,p)a i {k,p) + X)trlog[l + ^] +K 

7=1 k,p 1=1 



A contribution of the first multiplier in the integrand to the saddle point equations is of the 
factorizable form and it can be shown that this contribution can be neglected. So we can 
use again solution flTTp. Hence 



71-1 



(^(kMk)^} = \v(k)\ 2 \B\ 2 J2 fci+...+* n _ l +*,o II Hh)\ 2 = \v(k)\ : 

ki,...,kn-l 



i=l 



The last equality is due to the normalization condition (|3|). 
Finally we consider the average (|Hp. We have 

<^b| J24(k)4(-k + q)c 3 {-k' + qMk')\^ h } = (i/) h \c{(k)c{(-k + q)c n (-k' + g)c 1 (A; , )|^b) 



E 







d 



d 







VK*0i 2 EK?)f 



Tj=Tj=0 



, ; drj^k) dr) n (-k + q) drj n (-k' + q) drj^k' 
+C / DaDo- finiki, . . . , k n -i, -k + q)<p n (Pu ■ ■ -,Pn-u ~k' + q) JJ I 1 + ^il^ikhPi) 
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71-1 



fcl,...,fc n _l,g 
Pl.—.Pn-l 



7=1 



x[l + ai] x (A;i,A;)[l + ai] (pi, A;') exp 



n— 1 



n-l 



- EE^l^^^p) + E tr lo s[i + ^] + # 

7 = 1 fc,p 7=1 
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Again it can be shown as before that we can use saddle solution (|TT|), so we obtain 



J24(k)c}(-k + q) Cj (-k' + q) Ci {k')\^) = 6 k , k! \v(k)\ 2 T, K?)f 



n-1 



+v(k'HkHk')u(k')J2v(-k + q)v(-k' + q)\B\ 2 £ 5 k2+ ... +kn _ 1+q , Y[ \v(k t )\ 2 

q k2,...,k n -i i=2 

In the weak coupling limit \v(k)u(k)\ is non-zero only in a narrow region around the Fermi 
surface, so we can put \k\ = \k'\ = k F . Introducing the following notation 



f(k,k') = J2v(-k + q)v(-k' + q) 



|fc|=|fc'|=fc F 



we obtain (|14D. 
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